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Abstract. We describe our recent proposal of a path integral formulation of classical 
Hamiltonian dynamics. Which leads us here to a new attempt at hybrid dynamics, which 
concerns the direct coupling of classical and quantum mechanical degrees of freedom. This is of 
(■h ' practical as well as of foundational interest and no fully satisfactory solution of this problem has 

O |. been established to date. Related aspects will be observed in a general linear ensemble theory, 

which comprises classical and quantum dynamics in the form of Liouville and von Neumann 
equations, respectively, as special cases. Considering the simplest object characterized by a 
■ two-dimensional state-space, we illustrate how quantum mechanics is special in several respects 

3 ^ among possible linear generalizations. 
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' 1- The notorious linearity of the evolution equations 

This article is motivated by discussions of the origin of quantum mechanics, as presented, for 
example, in Refs. [U [2]. These suggest that quantum mechanics could possibly emerge by 
coarse-graining and sufficiently far from the Planck scale. Thus, the quantum phenomena would 
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arise from dynamics beneath. 

Numerous indications exist which support this view. Nevertheless, the general principles 
governing such hypothetical dynamics are under debate 01 El El El [10], [TH [121 [131 E] 
and, perhaps, have not even been touched upon, so far. 

In this situation, it may be worth while to focus on aspects of quantum and classical mechanics 
which bring them as closely together as possible, yet sharpen remaining differences. The linearity 



of the respective evolution equations must certainly be counted among such aspects. 

Indeed, it is remarkable that dynamical evolution in physics, to a large extent, can be encoded 
in one or the other of three linear equations: 

dtp = {H,p} , (1) 

id t ^ = , (2) 

idtp = [H,p] , (3) 

which, respectively, are the Liouville equation ([1]), governing the evolution of a classical phase 
space density p in terms of the Poisson bracket { , } and the Hamilton function H of the 
object under study, the Schrodinger equation ([2]) for the quantum mechanical wave function(al) 
in terms of the appropriate Hamilton operator H, and the von Neumann equation ([3|) for the 
density operator p in terms of the commutator [ , ] with H. 



We emphasize that all three equations pertain to (or can be interpreted as) ensemble theories. 
They have the generic structure: 

id t "state" = t "state" , (4) 

where the "state" appears only linearly in all cases. This first-order in time differential equation, 
of course, assumes widely different forms depending on the forces or dynamics under study. That 
is, the "Liouville superoperator" C on the right-hand side allows to incorporate physics ranging 
from classical non-relativistic single particles all the way to relativistic quantum fields, employing 
a phase space ensemble of initial conditions in classical statistical mechanics and the functional 
Schrodinger picture for interacting quantum fields. 

Since the Schrodinger equation describes special cases (pure states) of situations (pure and 
mixed states) covered by the von Neumann equation, we consider Eqs. ([1]) and ([3]), eventually in 
the form of Eq. Q, in the following. 

In passing, we recall that nonlinear modifications of quantum mechanics have been proposed 
and severely constrained in various ways, see, for example, Refs. [T5 l [TB I HT j US]. 

In distinction, modifications of the linear Eq. ([3]), in particular, have been derived for "open 
systems" or quantum mechanical objects interacting with an environment (such as a heat bath 
of oscillators, etc.). - Here, one is clearly not interested in modifying quantum mechanics but 
in applying it to complex situations. - The generic form of the resulting consistent dynamics, 
replacing Eq. ([3]) while preserving the defining properties of the density operator, is known as the 
Lindblad equation [19, 20j. Its generator provides an important example of a completely positive 
map. Some aspects of this will be discussed in more detail in Section 4, in order to contrast it 
with modifications of the linear dynamics introduced there. 

A hypothetical fundamental coupling of quantum mechanical and classical degrees of freedom 
- "hybrid dynamics" - instead, presents a departure from quantum mechanics. While this does 
not touch the linearity, in the sense of Eq. 0, consistency of the resulting theory has to be 
carefully examined. 

We remark that there is also a practical interest in certain forms of hybrid dynamics. The 
Born-Oppenheimer approximation, for example, is based on a separation of interacting slow and 
fast degrees of freedom of a compound object. The former are treated as approximately classical 
while the latter as of quantum mechanical nature. It must be emphasized, however, that in 
these cases, hybrid dynamics is considered as relevant approximate description of an intrinsically 
quantum mechanical object. This has recently been employed, for example, in a re-derivation 
of geometric forces and Berry's phase [22] . 

Such considerations are and will become increasingly important for precise manipulations of 
quantum mechanical objects by apparently, for all practical purposes classical means, especially 
in (sub-)nanotechnological devices. 

However, this may be also relevant for studies of the measurement problem: Namely, quantum 
theory, endowed with the Copenhagen interpretation, assumes a coupling of the quantum 
mechanical object to the classical measuring (or manipulating) apparatus. However, it remains 
profoundly silent about how this is realized and how the ensuing state reduction ( "collapse of the 
wave function") proceeds. - Truly hybrid dynamics, which does not arise from an approximation, 
augmented by effects of environment induced decoherence and anharmonic forces, must have 
something to say about measurement situations, if it can be formulated consistently. 

Presently, we are motivated by interest in the back-reaction effect of quantum fluctuations on 
classical degrees of freedom, in particular if they are physically distinct. We recall here especially 

1 Mean field theory, based on the expansion of quantum mechanical variables into a classical part plus quantum 
fluctuations, leads to another approximation scheme and another form of hybrid dynamics. This is discussed more 
generally for macroscopic quantum phenomena in Ref. [21| . 



discussions of the "semiclassical" Einstein equation coupling the classical metric of spacetime 
to the expectation value of the energy-momentum tensor of quantized matter fields. Can this 
be made into a consistent hybrid theory leaving gravity unquantized? This has recently been 
re-examined in Ref. [23]; earlier related work includes Refs. [Ml 1251 1261 127112811291 1301 IHT] . 

Concerning the origin of quantum mechanics from a coarse-grained deterministic dynamics 
[1, 2j, to which we alluded above, the back-reaction problem can be more provocatively stated 
as the problem of the interplay of fluctuations among underlying deterministic and emergent 
quantum mechanical degrees of freedom. Or, in short: "Can quantum mechanics be seeded?" 

In this paper, we present a few ingredients which may be helpful in further studies of these 
questions. - In the following Section 2, we begin with a brief recapitulation of our path integral 
formulation of classical Hamiltonian dynamics. Some simple observations following from this 
lead us to take up the subject of hybrid dynamics in Section 3. In Section 4, we continue the 
discussion from a different angle, namely with some remarks on problems encountered when 
trying to embed classical or quantum mechanics in a more general linear dynamics. 

2. A path integral for classical Hamiltonian dynamics 

Our aim here is to solve the Liouville equation ([1]) with the help of a suitable propagator which, 
in turn, is represented as a path integral. We summarize essential steps, while more details of 
the derivation can be found in Ref. [32] , 

The following considerations are rather independent of the number of degrees of freedom and 
apply to matrix or Grassmann valued variables as well; field theories require a classical functional 
formalism, which has been considered elsewhere [HI [33]. We consider a one-dimensional system, 
for simplicity. 

2.1. The quantum-like version of the Liouville equation 

We assume conservative forces acting on the classical object under study and that Hamilton's 
equations are determined by the generic Hamiltonian function: 

H(x,p):=^p 2 + V(x) , (5) 

where x and p denote generalized coordinate and momentum, respectively, and where V(x) 
stands for an external potential (a mass parameter will be inserted later). 

An ensemble of such objects can be described by a probability distribution function p 
depending on the x, p-coordinates of phase space and time. Such a distribution evolves according 
to the Liouville equation: 

-**=£-!-£-g=w.-™»,. w 

with V'(x) := dV(x)/dx. - The relative minus sign in the Poisson bracket, or between terms 
here, reflects the symplectic phase space symmetry. It will give rise to a commutator structure, 
which reminds one of quantum mechanics, as we shall see momentarily. 

Following a Fourier transformation, p(x,p;t) = f dy e~ tpy p(x,y;t), and a transformation of 
the effectively doubled number of spacelike coordinates, 

Q:=x + y/2 , q := x - y/2 , (7) 

the Liouville equation becomes: 



idtp = {H Q -H q +£(Q,q)}p , (8) 



H x := -\d* + V{x) , for X = Q,q , (9) 

£(Q,q) := (Q-q)V'(^-)-V(Q) + V(q) = -£(q,Q) . (10) 

Thus, it bears strong resemblance to the von Neumann equation, cf. Eq. ([3]), considering g; t) 
as matrix elements of a density operator />(t). 

We automatically recover the Hamiltonian operator H related to the Hamiltonian function, 
Eq. dSJ), as in quantum theory. Furthermore, reality and normalization of the phase space 
probability distribution p(x,p; t) translate into hermiticity and trace normalization of the density 
operator p(t) [32J. 

However, an essential dynamical feature here consists in the interaction £ between bra- and 
fcei-states. Thus, generally, the Hilbert space and its dual are coupled by a genuine superoperator, 
a concept to be defined in the following subsection. 

It is remarkable that the interaction between bra- and ket-states vanishes under certain 
circumstances: 

£ = <^=^ potential V(x) is constant, linear, or harmonic , (11) 

rendering a Liouville superoperator of quantum mechanical form, i.e., as in Eq. ©.0 This has 
been discussed in Ref. [2] under the perspective of having quantum phenomena emerge due to 
discrete spacetime structure. 

In the following, we will study in more detail the classical dynamics described by Eq. (JSD, or 
by appropriate generalizations, and pay particular attention to the presence of the superoperator 
£ , when comparing with the von Neumann equation and its solution by a propagator. 

2.2. The super space 

The dynamics of density operators in the general form of Eq. Q can be conveniently rewritten 
by introducing the concept of superspace, also called Liouville space [32| 134"]. Fl 

Let H denote a Hamiltonian operator, as in quantum theory, and let {\j)}, j = 1, . . . , N, 
present a complete orthonormal set of basis states, assuming here that the relevant Hilbert space 
is iV-dimensional. Then, the matrix elements of the von Neumann equation (J3j) read: 

idtp jk = [{Hp) jk - {pH) jk } , j, k = 1, 2, . . . , N , (12) 
with a density matrix p of iV 2 elements. Or, written as in Eq. @: 

idtPjk = 22 £jk,lmPlm = (£-P)jk ) 
l,m 

where the Liouville superoperator C is now defined by its matrix elements: 

£jk,lm '■= HjlO~mk — SjlH m k . (14) 

2 The interaction £ is antisymmetric under Q o q. It follows that the complete Liouville superoperator on the 
right-hand side of Eq. (J8J, to be compared with Eq. ((4]), has a symmetric spectrum with respect to zero and, in 
general, will not be bounded below. Related observations were discussed, for example, in Refs. [51 [51 1^1 ITU]. 

3 Analogously, the vanishing of £ in a field theory amounts to having massive or massless free fields, with or 
without external sources, and with or without bilinear couplings. In these cases, anharmonic forces or interactions 
are absent. 

4 No relation with the case of supersymmetry is implied. 



(13) 



This simple rewriting may suggest to introduce a space in which the density operator is a vector. 
This is the role of the Liouville space (or superspace). 

The dynamics of density operators can then be described in parallel for classical and quantum 
mechanics. In general, they will differ, of course, by the precise form of the superoperator, 
examples of which we have seen here and in the previous subsection. 

Given the Hilbert space, as above, the density operator can be expanded as: 

p = T,pokb)^\ ■ ( 15 ) 

The family of N 2 operators \j)(k\, with j, k = 1, ■ ■ • , N, can be interpreted as a complete set of 
matrices, or vectors, such that the density operator becomes: 

\p) = T,pjk\m , (i6) 

3,k 

where the "ket" \jk)) denotes the Liouville space vector representing the Hilbert space operator 
\j){k\. Similarly, we introduce a "bra"-vector {{jk\ as the Hermitian conjugate of \jk)). 

Consequently, any operator A is represented by a vector, denoted by \A)), and can be 
expanded as: 

\A))=Y J A ]k \jk)) , (17) 

j,k 

where Aj k are the usual matrix elements {j\A\k). - Furthermore, with a bra- vector ((B\ 
representing B\ the scalar product of two operators is defined by: 

((B\A)) := Tr(B^A) . (18) 

This implies the orthonormality condition: 

({jk\mn}} = Tr(\k)(j\m){n\) = S kn 8 jm , (19) 

in analogy to and based on (j\k) = 5jk- - Finally, consider the scalar product: 

((jk\A)) = Tr(\k)(j\A) = -£(l\k)(j\A\l) = (j\A\k) ee A jk . (20) 

I 

Upon substitution in Eq. (fTT|) . this yields: 

\A)) = J2\jk))((jk\A}) . (21) 

3,k 

This is consistent with the following completeness relation in Liouville space: 

J2\jk))((jk\ = l . (22) 

3,k 

It is easy to see now that the Liouville space is a linear space and that the density operator p, 
in particular, is a vector in this space. - To conclude these formal considerations, we define a 
linear operator in terms of its matrix elements: 

F ■■= \jk)){{jk\T\mn)}{(mn\= ^ ^jk,mn\jk))((mn\ , (23) 

j,k,m,n j,k,m,n 



an example of which is given in Eq. (|14p . 

The importance of Liouville space for classical and quantum dynamics resides in that Liouville 
and von Neumann equations, both, assume the form of Eq. incorporating an appropriate 
superoperator C, cf. Eqs, ([8])-(fT0]) and (fT2]) - ([T^|) , respectively. Thus, we find a close formal 
similarity between the structure of these equations and the Schrodinger equation ([2]). 

Therefore, it is plausible that formal derivations or techniques concerning the solution of the 
Schrodinger equation can be transferred to the cases of Liouville or von Neumann equations with 
the help of Liouville space notions. This regards perturbation theory as much as nonperturbative 
methods, the path integral approach in particular, to which we turn next. 

2.3. The Liouville propagator and path integral 

The technical ingredients of the Feynman path integral approach, for the derivation of quantum 
mechanical propagators in particular, are very well known |35l 136] . We will make use of these 
ideas, in order to derive a path integral for the propagator of density matrices based on the 
Liouville space formulation of the preceding subsection. 

Our derivation relies on the close formal similarity between the classical Liouville equation 
and the von Neumann equation on one hand side and the Schrodinger equation on the other, in 
the appropriate representation that we discussed. El 

In particular, the formal solutions of the classical Liouville equation and of the quantum 
mechanical von Neumann equation, both, can be written in the form: 

\p(t))) = e -^C*-*o)/R| p(to)>> , (24 ) 

where C is the relevant Liouville superoperator. Here, we have: 

«Q, q\£\Q', q')) = 5(Q - Q')5(q - q')(H(Q) - H(q) + £(Q, q)) , (25) 

where H denotes the appropriate Hamilton operator in coordinate representation, as indicated, 
which alone is relevant for the von Neumann equation, while £ represents the additionally present 
superoperator term for classical dynamics, cf. Subsection 2.1. 

In order to solve the problem of time evolution in the present case, we need to know the 
(super)matrix elements entering the propagation equation: 

((Q,q\p(t))) = J dQ'dq' ((QM^ it -* )/h \Q'M((Q',J\pto))) , (26) 

in analogy to the case of a state vector evolving according the Schrodinger equation. Not 
surprisingly, we can now follow the usual steps [35} [36] , in order to construct the path integral 
representation of our propagator [32]. 

In this derivation, one has to pay attention to a suitable generalization of the Trotter product 
formula. This works out in a straightforward way; the relevant definitions and details of the 
proof are given in Ref . |37] . 

Rewriting the Eq. ([26]) as: 

p(Q, q;t) = J dQ'dq' g(Q, q - t\Q', q '- t ) P (Q', q '- 1 ) , (27) 

our interest is to know the superpropagator Q. 



5 In this subsection, we reinstate h explicitly. 



Following the recipe to arrive at a path integral representation, we make use here of suitably 
inserted complete sets of superspace vectors, such as: 

J dQdq\Q,q))((Q,q\ = l , (28) 

and, correspondingly, for momentum space, cf. Eq. (|22p . Using the plane wave relation between 
coordinate and momentum eigenfunctions, we also employ: 

((P,p\=^ K fdQdqexp(-^(PQ-pq))((Q,q\ . (29) 

Furthermore, the orthogonality relation {{Q,q\Q' ,q')) = 8(Q — Q')5(q — q'), cf. Eq. (fT9|) . implies: 

((P,P\Q, q)) = ^ exp ( - % - (PQ - pq)) . (30) 

Then, with all following steps of the derivation in parallel with the usual ones, it is 
straightforward to obtain the Liouville path integral [32, 37J: 

G(Qf,qf,t\Qi,Qf,to) = j VQVq exp(~S[Q,Q;q,q]) , (31) 

with the boundary conditions Q(U) = Qi, q(U) = qi, Q(tf) = Qf, and q(tf) = (?/, otherwise 
unrestricted paths, and with the superaction S defined by: 

S := f dr (f Q 2 - V(Q) - (f q 2 - V{q)) - £(Q, q)) , (32) 
Jto 

for a particle of mass m. We recall that £ = corresponds to evolution according to 
the von Neumann equation, whereas £ ^ represents the only modification due to classical 
dynamics, in accordance with the Liouville equation, cf. Eqs. ([8|)-(fT0]). 

However simple this result may seem, the Eqs. (|3ip - (|32p describe time evolution of the full 
density matrix. The particular new feature is that formally classical dynamics is treated on 
the same footing as quantum mechanics, differing only in the action entering the phase in the 
integrand of the path integral. 

These considerations are limited neither by one dimension nor by single-particle physics, but 
can be extended all the way to relativistic field theories. This offers new calculational tools, 
new approximation methods in particular, and may be of interest for applications in classical 
statistical mechanics. In the following, however, we turn to the quantum-classical divide. 

2.4- Simple properties of the Liouville path integral 

First of all, it seems worth while to record a few simple properties of the Liouvile path integral. 

We observe that in the absence of forces, V = 0, the classical and quantum mechanical 
propagators coincide, according to Eqs. (f3Tj) - ([32"j) . Which implies that classical and quantum 
mechanical behaviour can differ at most in the initial states that are being propagated, in this 
case. - This holds true even for V ^ 0, if £ = 0, i.e., for harmonic forces, cf. (jlip . - Thus, 
the cherished "textbook effect" of the spreading of wave packets is not a peculiar (kinematical) 
quantum effect, but rather an effect of the particular states considered! 

A simple calculation, taking into account our above transformations from x,p- to Q,q- 
coordinates (and back), shows that a classical point particle with initial phase space distribution 
oc S(x — xo)5(p — pq) is propagated to the distribution oc 5{x — xq — pot/m))5(p — po), i-e. along 



a straight line path, as expected. - Instead, a free Gaussian wave packet in Q, g-space spreads 
in the way described in textbooks on quantum mechanics. 

Finally, we can convince ourselves that the classical dynamics is properly represented by the 
Liouville propagator, in general, by undoing the coordinate transformations ([7]) directly in the 
path integral representation of the superpropagator: 

y:=Q-q , x := ^y-^ , (33) 

a transformation with unit Jacobian. Following a partial integration of the superaction, the 
path integral over the y-coordinate simply yields a functional 5- "function" and results in: 

Q(xf,yf,t\xi,yi]to) = J Vx exp(imxfyf) 5[mx + V'(x)} ex.p(-imxiyi) , (34) 

with the boundary conditions x(tij) = Xij, similarly for the velocities, and where yij denote 
the initial and final values of the y-variable, respectively. The phase factors stem from the partial 
integration of the superaction. 

Thus, we find the expected result that only paths following solutions of the classical equations 
of motion contribute to the propagator. The role of the phase factors is easily understood by 
recalling Eq. (J27J): After the inverse coordinate transformations ([33]), one of the two ordinary 
integrations there becomes here, including the relevant (initial) phase factor: 

J dyi e~ tmXlVl p(xi, yi\ t ) = p(xi,mxi = p»; to) > (35) 

i.e., this incorporates the inverse Fourier transformation, back to the momentum variable of 
phase space. The second (final) phase factor, then, is necessary for the propagator to fulfil the 
important semi-group property [35J. 

3. Hybrid dynamics for two interacting objects 

We have seen that the path integral for the propagator of the classical Liouville equation shows 
significant similarity with the propagator of the quantum mechanical von Neumann equation. 
This suggests a new perspective on hybrid dynamics - the hypothetical direct coupling of 
quantum and classical degrees of freedom that we discussed in Section 1. 

In particular, we are interested here in the propagator for the density matrix of a bi-partite 
system, composed of a classical and a quantum mechanical particle of masses m, m! in external 
potentials V, v, respectively. Correspondingly, the relevant coordinates will be denoted by Q, q 
and Q',q', respectively. Our previous considerations lead us to propose the following hybrid 
superpropagator (cf. Eq. (|3T])): 

g h (Qf,q f ;Q' f ,q' f ;t\Q i ,q i ;Q' i ,q' i ;t )= J VQVqVQ'Vq' exp (^S[Q,Q; q,q;Q' ,Q'; q' ,q'}) , 

(36) 

where the superaction naturally consists of three contributions, S = 5 c j + 5q m + 5^: 



S d := £dr(f(Q 2 -g 2 )-(Q- (7 )y'(^)) , (37) 

%n := f dr (f (Q> 2 - q' 2 ) - (v(Q') - v(q'))) , (38) 

S h := - f* dTV(Q,q;Q',q') . (39) 
J to 



Several remarks are in order here. The contribution 5 c i is the classical superaction arrived at 
in Eqs. (j31~]) - (f32j) : we just inserted £ , Eq. (fl~0|) . explicitly and collected terms. Correspondingly, 
the contribution 5qm is the action for the propagator of the von Neumann equation, i.e., it 
describes the quantum mechanical particle in the usual way. Thus, in the absence of the 
interaction term 5^, the hybrid superpropagator Qy. factorizes into the corresponding ones for 
Liouville and von Neumann equations. This would consistently describe a composite system of 
two independent, classical and quantum mechanical particles. 

The interaction term introduces the coupling responsible for hybrid dynamics. - For our 
present purposes, we restrict this coupling by two consistency requirements: |j 

(A) 5^ + 5^ describes a classical particle interacting with another particle (coordinates Q',q'), 
as if the latter were classical; 

(B) Sqm + 5^ describes a quantum mechanical particle interacting with another particle 
(coordinates Q,q), as if the latter were quantum mechanical. 

In view of the general form of the external potential terms in Eqs. (|37p and (|38p . these 
conditions can only be fulfilled, if V(Q,q;Q',q') is harmonic, i.e., is a polynomial of degree 
less than or equal to two in all variables. Thus, the hybrid coupling is of a form that admits 
a classical or quantum mechanical interpretation, depending on whether it is viewed from the 
classical or quantum mechanical subsystem. 

An example is provided by a distance dependent oscillator potential, A _1 V := (x — x') 2 . 
Following the derivation of the classical superpropagator, this becomes, in the above coordinates: 

A- 1 V(Q, q; Q',q') = ((Q - q)d {Q+q)/2 + (Q' - q')d {Q , +q , )/2 ) ((Q + q)/2 - (Q' + q')/2) 2 

= (Q-q)(Q + q) + (Q'-q')(Q' + q') 2 -2(QQ'-qq / ) , (40) 

i.e., it can be seen as coupling between two classical particles, fulfilling condition (A). On 
the other hand, the separable oscillator terms here also equal Q 2 — q 2 + Q' 2 — q' 2 , i.e., are of 
quantum mechanical form, in accordance with (jlip . and the bilinear coupling can as well be 
seen as quantum mechanical, fulfilling condition (B). 

Interesting consequences of these conditions, determining a harmonic hybrid interaction, are 
summarized in the following Table: 
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Indicated are the nature of the potentials acting on the classical and quantum mechanical particles and 
their interaction, V, v, and V, respectively - h: harmonic; anh: anharmonic - and the character of the 
resulting hybrid dynamics - CL: classical; QM: quantum mechanical. 

For example, if the classical and quantum mechanical particles are governed by a harmonic 
potential and an anharmonic potential, respectively, then the resulting dynamics of the 
composite system can be considered as that of two bilinearly coupled quantum mechanical 
particles, one of which moves in a harmonic potential. - The case of both potentials being 
anharmonic is most interesting. However, no general statements can be made in this case 
without further study. The quantum mechanical or classical behaviour of the composite system 
might very well depend on where (concerning its variables) one is looking at this object! Lj 

6 The question of consistency of the map induced by the superpropagator which evolves the quantum-classical 
hybrid density will be discussed in more detail elsewhere. 

7 These findings might help with the problem posed in Section 1: "Can quantum mechanics be seeded?" 



3.1. Intra- and interspace entanglement 

There is a qualitative difference between CL and QM, contained in the full path integral for the 
superpropagator, Eqs. ()36p -(|39 p . Following standard arguments which lead back from a path 
integral to an equivalent equation of motion for the (hybrid) density operator p, several related 
observations may be interesting. 

The QM evolution is generated by a commutator of the Hamiltonian with the density 
operator. Generally, this superposes and, for multi-partite systems, in particular, entangles 
underlying bra- and ket-states separately, oc HijPjk — PijHj^ (using a convenient notation with 
discrete indices). For a bi-partite system, it is revealing to write the relevant interaction terms 
explicitly: 

[Hint, p] = (8) H 2 p 2 - piHt (8) p 2 H 2 , (41) 

for an interaction Hi n t := H\ ® H 2 , with the factors acting on subsystems "1" and "2", 
respectively, and where p = p\®p 2 , for a separable initial state. This has been called dynamically 
assisted entanglement generation, see, for example, Refs. [38j [39j 00] . 

It may come as a surprise that the CL evolution does this just as well, due to the structure of 
the superoperator. For polynomial interactions, in particular, the superoperator always contains 
a contribution proportional to the usual QM terms. 

However, the CL evolution, generally, produces additional correlations in p, due to terms 
contained in Cij-tiPki which entangle bra- and ket-states. - In comparison with Eq. ()4ip . such 
terms can have the unfamiliar structure: 

H'^® p 2 H' 2 - p!^(8) H' 2 p 2 , (42) 

which differs decidedly from a commutator. - This leads us to distinguish intra- (i.e., within 
given tensor product Hilbert space of subsystems "1" and "2") and inter-space entanglement 
(i.e., between said Hilbert space and its dual). 

Consider, for example, the anharmonic potential V(x\ — x 2 ) := A(xi — x 2 ) 4 for a bi-partite 
system consisting of particles "1" and "2" . Similarly as before, this leads here to the interaction: 

V(Qi, qr, Q 2 , g 2 ) = \\{Qi -qi- (Q2 - <a))(Qi + gi - (Q 2 + ^)) 3 , (43) 

in terms of, by now, familiar variables, taking into account subsystems "1" and "2" with Q's 
and q's refering to bra- and ket-states, respectively. Besides separable terms, ex (Q a — q a ){Qa + 
g a ) 3 , a = 1,2, there are terms which mix (and entangle) variables of both subsystems, as 
usual in QM. However, there are also additional terms that refer to Hilbert space and its dual 
simultaneously (and entangle corresponding states), for example, cx QaQbqti b ^ a. 

In retrospect, somehow, such difference between CL and QM evolution had to be expected: 
instead with superstates \Q,q)), we could have worked with superstates \x,p)), relating to 
coordinates and momenta of the classical theory. There, coordinates and momenta end up 
tightly correlated, due to Hamilton's equations, and produce inter-space entanglement in an 
interacting bi-partite system. 

Thus, the confrontation of CL with QM, as in our side-by-side study, is quite revealing. In 
particular, we speculate that this opens new views on generating entanglement in multipartite 
systems, perhaps, by evolving through quasiclassical stages or by making use of decohered 
intermediary states. |f| 

8 Previous considerations of the semiclassical regime, such as in Refs. [381 139] , were motivated as suitable 
approximations of the quantum mechanical evolution, in particular, for studies of the different decoherence 
properties between classically regular and chaotic systems. Our results seem to show that crossing the quantum- 
classical divide may offer an additional resource for entanglement generation and related "truly quantum" 
phenomena. This might be related to a common "underlying reality" of CL and QM physics, assumed to consist, 
for example, only in statistical correlations in Refs. [411 142] . 



Concerning the quantum-classical divide, the present analysis shows that there is an 
appealing, if not puzzling formal similarity between CL and QM. However, this demonstrates 
one more time that what has been discussed in various ways as CL limit of QM - and which is 
similarly relevant for "emergent QM" - deserves further study. 

While our work has been concerned mainly with the evolution of CL or QM objects, we recall 
that V.I. Man'ko and collaborators have pointed out that classical states may differ widely from 
what could be obtained as the a h — > 0" limit of quantum mechanical ones. They show that all 
states can be classified by their 'tomograms' as either CL or QM, CL and QM, and neither CL 
nor QM 03]. 

The classical limit can be considered a limit "For^UlPracticalPurposes" , gradually approached 
with decoherence as an essential but insufficient ingredient or, formally, following the mnemonic 
u h —> 0" rule. However, in order to truly bridge the qualitative difference between intra- and 
inter-space entanglement that we find, and explain the "Man'ko classes of states" , some unknown 
dynamics seems missing, d 

The problems discussed here lead us to the question: "Does the h — > deformation of QM 
provide the only interesting linear dynamics besides QM itself?" 

4. How special is quantum mechanics? 

We are interested here once more in the structure of the linear evolution equations that we have 
discussed side by side, namely the classical Liouville equation and the von Neumann equation 
of quantum mechanics. However, we consider this with respect to a hypothetical most general 
linear dynamics that can be represented in the generic form of Eq. (j5|): 

9 tPij = ~Y1 ZijMPkl , ( 44 ) 
1 k,l 

where the indices k, I = 1, . . . , N) refer to a discrete, finite dimensional Hilbert space, which 
we assume for simplicity. Thus, generally, there are iV 2 x N 2 complex coefficients to be specified. 
Imposing the constraint that the 'density matrix' is and remains Hermitian, which requires that 
£-ij,M = — ik, the set of coefficients can be specified by iV 4 real parameters. 

This should be compared with the von Neumann equation - which we discussed, in order 
to introduce the concept of a superoperator, cf. Eqs. ([T2|) - ([TS]) in Section 2.2. In terms of a 
Hermitian Hamiltonian, H, the matrix elements of which are specified by N 2 real parameters, 
the corresponding superoperator has been obtained as: Aj.iW = Hik^ij — SikHlj- 

Therefore, we can write a symbolic relation, concerning the number of real parameters which 
determine the respective dynamical equation: 

"QM~\/GL" , (45) 

i.e., the number of real parameters entering the quantum mechanical evolution law (QM) scales 
with the Hilbert space dimension like the square root of the corresponding number for the most 
general linear dynamics (GL) which preserves hermiticity. By choosing the eigenstates of the 
Hamiltonian as basis of the Hilbert space, the number of relevant parameters in QM could be 
further reduced to N, the number of real eigenvalues of H. 

We may wonder whether the commutator in Eq. (|I2p . [H, .], presents the minimal structure 
preserving hermiticity, normalization, and positivity of the density matrix. - Conversely, can 
the different numbers of real parameters between a more general linear evolution law and QM 
be attributed to an attractor mechanism or information loss if and when QM is emergent [HE]? 

9 A simple attractor model, motivated by assumptions about effects of fundamental spacetime discreteness [2], 
has been discussed in Ref. [3]. 



As a first step to investigate these issues, we consider here whether QM admits a more general 
linear evolution. - We recall the following result for the QM of open systems, see Ref. [H] and 
references there: 

• If and only if the Liouville superoperator C entering the right-hand side of Eq. (|44p generates 
a Hermitian, trace preserving, completely positive map, M.(t)p(0) := exp(— iCt)p(0), then 
it can be written in canonical Kraus form: 

M(t)p(0) =Y,Mk(t)p(0)M k \t) 

k 

with EfcM^Mfe = l NxN . 

• Such a map Xi applied to a Hermitian, positive semidefinite, normalized density operator 
yields another one and linear maps beween density operators can be written in Kraus form. 

An evolution equation corresponding to such a map, in general, contains the von Neumann 
equation as a special case. - In fact, the solution of the von Neumann equation (|12p is provided 
by a unitary transformation, p(t) = U (t) p(0)U (ty , which presents the simplest case with only 
one unitary Kraus operator, M(t) = U(t) := exp(-iHt). 

Therefore, any generalization which maintains the defining properties of density operators, 
generally, will need additional Kraus operators specified by additonal parameters. - In this sense, 
the commutator defining the von Neumann equation presents the minimal structure indeed. 

Furthermore, it has been shown that the dynamics generated by M of the Kraus form, in 
general, is equivalent to a non-unitary reduced dynamics of a unitary dynamics on a bigger 
(tensor product) Hilbert space "System (g> -Environment" [44] : 

Mp s = Tr^ (U b igps pEUlg) 

with UbigUlig = ul ig Ubig = 1. Hence, we learn: 

• Hermitian, trace preserving, completely positive maps, while generalizing the von Neumann 
dynamics, are not general enough to "leave QM". They describe the QM of open systems. 

However, this kind of maps does not exhaust the larger set of maps generated by all possible 
Liouville superoperators. 

We conclude that one cannot "leave QM" by invoking more general linear dynamics than 
that generated by trace preserving, completely positive maps without affecting properties and 
interpretation of the states represented by density matrices. 

It is useful to recall here (the motivation behind) the assumptions made concerning the 
properties and interpretation of density operators. - A density operator p is required to be i) 
Hermitian, i.e., to have real eigenvalues, in order to qualify as an observable. This is needed, in 
turn, if one requires p to be ii) positive- semidefinite and iii) normalized. All three properties, 
together, are necessary assumptions for the standard probability interpretation of the eigenvalues 
of a density matrix, according to the Born rule. 

A recurrent theme, when comparing quantum with classical states, for example, with the 
help of the Wigner (function) transform of the density matrix, is the appearance of negative 
eigenvalues or "negative probabilities". Numerous attempts have been made to give a satisfactory 
physical interpretation and mathematically consistent definition to these, see Refs. [42|. |4"5| |4"6] 
and further references there. 

Instead of entering this discussion, we look at the arguably simplest example of general linear 
dynamics in the following. This is obtained by giving up the requirement ii) above, i.e., by 
abandoning positivity, and by considering a two-dimensional Hilbert space. 

10 It will be interesting to see how (any form of) hybrid dynamics, cf. Section 3., fares in this respect. 



(46) 



(47) 



4-1- General linear 2x2 dynamics: a model 

A convenient parametrization of the most general Liouville superoperator C for a two- 
dimensional state space can be written with the help of the Pauli matrices a. Thus, the 
(super)matrix elements Cijkh i,j,k,l = 1,2, are defined by: 

Cij,kr-=iG»AA a » + b »-v\ k )( l \ a '' + K-0\j} , (48) 

where a summation over repeated indices is understood; in this notation we have (Cp)ij = 
£>ij,klPkl- Hermiticity requires = G*^. Generally, three vectors 6^=0,1,2 £ R- 3 ar e needed. 
Furthermore, we define: 

/I \ 

Q := 011 5-12 (49) 
V g* 2 922 J 

i.e., where g is a Hermitian 2x2 matrix. This Ansatz saturates the expected number of 
iV 4 = 2 4 = 16 real parameters, since a M , 6 M , and g contribute 3, 9, and 4 parameters, respectively. 

In order that the map generated by C be trace preserving, we must have = idtTvp = TrCp, 

for all p, which leads to Cu^i = 0. This yields four real constraints, 

a-Q -a + b-G -b = 0, (50) 

2a ■ Reg ■ b + igijbi xbj = , (51) 

which reduce the number of available real parameters to twelve. 

Next, we conveniently choose the vectors 60,1,2 to form a right-handed orthogonal system; 
with respect to suitable coordinates, this sets six vector components to zero. Thus, we are left 
with six real parameters, three of which could characterize the QM evolution of a two-state 
object, while three pertain to the generalization we are concerned with. 

Finally, we assume for simplicity that the 'QM part' of C is diagonal, corresponding to a spin- 
1/2 particle in a constant external magnetic field parallel to the quantization axis, for example. 
Setting one of the remaining parameters to zero, we obtain a simplified model with altogether 
three real parameters, a, (3, 7 G R: 

Aj,fei : = a^(i\a z \k)5ij - 5 ik (l\a z \j)^ 

+Pl{{iWz\k)(l\o- x \j) - (i\a x \k)(l\<T z \j)} + i^^ay^Sij + 5 ik (l\a y \j}} 

+if3 2 5 lk 5 l3 +i 1 2 ii\a y \k){l\a y \j) - i((3 2 + 7 2 )(i\a z \k)(l\a z \j) , (52) 

where QM terms oc ao~ z feature in the first line. All other contributions on the right-hand side 
have no counterpart in QM; this will become obvious by explicitly solving the model. The first 
two terms in the second line and the two last ones in the third couple bra- and ket-states, which 
we discussed as a consequence of classical evolution in Sections 2.1 and 3.1. It is easy to verify 
that the generator — iC is Hermitian and trace preserving, as it should. 

Taking hermiticity and trace normalization of the density matrix into account by: 

p=( P f / 12 ) , (53) 
V P12 1 - Pn J 

we write the resulting evolution equation explicitly: 

id tP = Cp = (54) 

i 7 2 (l - 2 Pu ) 2/? 7 + (2a + i[2(3 2 + 7 2 ])pi 2 - ij 2 p* 12 

-207 - (2a - i[2(3 2 + 7 2 ])pt 2 - i^ 2 p l2 -i 7 2 (l - 2p u ) 



Note that for /J, 7 —> 0, we recover the von Neumann equation, e.g., for a spin-1/2 particle 
in a constant magnetic field. In this case, besides the Hamiltonian, H := aa z , also the total 
spin-squared S 2 := <r 2 /4 oc 1 is a constant of motion. 

4-2. Constants of motion 

We remark that if a constant operator, C, obeys Tr(C/S(0)) = Tr(Cp(i)), for all solutions p of 
the general linear evolution equation (|44p . then CijCji ki = 0, which generalizes the vanishing of 
the commutator [C, H]ik in QM, for constants of motion. This is fulfilled for any C = cl, with 
constant c, due to the preservation of the trace normalization of the density matrix, incorporated 
by Cu m = 0- However, this raises also the interesting question, whether there exist conserved 
superoperators, Q, which can be defined by the vanishing supercommutator: 

\Qi£\ij,mn := Qij,kl£-M,mn £-ij,klQkl,mn = • (55) 

In the present case, this amounts to a 4 x 4 matrix equation, to be studied. 
4-3. Solution of the 2x2 model 

We now turn to the explicit solution of Eq. (|54|) . which can be represented as follows. The matrix 
elements are: 

Pn(t) = 1 - p 22 (t) = l + ( pil (0)-i)e- 2 ^ , (56) 
Pl2(t)=P2i(t) = r + e Q +* + r_e n "* + r c , (57) 
where r± are determined by initial conditions. The remaining constants are defined by: 

a - i(3 2 

rc := ^V+W + 7 2 ) ' (58) 
Q± := 2/3 2 + 7 2 ±2«a(l-7 4 /4a 2 ) 1/2 , (59) 

and we may choose a > 0. Diagonalizing the resulting traceless Hermitian density matrix, we 
obtain the eigenvalues: 

P± (t) = l± (( Pll (0) - l/2) 2 e- 4 ^ + \p l2 (t)\ 2 ) 112 , (60) 



with \p 12 \ 2 given by Eqs. (f5Tj)Hf59]) . 

Several features of the solutions of our model seem unexpected. - To begin with, for /3, 7 — )■ 0, 
we obtain: 

Pn(t) = 1 - p 22 (t) = Pll (0) , (61) 
Pi2(t) = p* 2l (t) = r + e 2M + r.e- 2iat , (62) 

which amounts to the usual QM solution only, if we choose r + = 0. From a QM perspective, 
the larger set of solutions, parametrized by a larger set of initial conditions, is quite surprising. 



11 Interestingly, the underlying equation for the off-diagonal matrix elements, for 7 7^ in particular, is of second 
order; it decouples into two first order equations in the limit 7 — > 0, independently of the value of {5. 



It leads to time dependent 'probabilities', if we try to interpret the eigenvalues of p in the usual 
way: 

p ± (t) = \ ± (( Pll (0) - 1/2) 2 + |r + | 2 + |r_| 2 + 2Re(r + r!.e 4fa <)) ^ . (63) 

While these 'probabilities' are real and normalized, by construction of our model, they might 
temporarily fall outside the interval [0, 1], depending on the parameters r±. 

This remark applies also for the general case, with j3, 7 7^ 0, where amplitudes grow 
oc exp(2/3 2 + 7 2 )i, see Eqs. (^7D .(l59l). 

The effect of such anomalous 'probabilities' is clearly seen in the following mean square 
deviations sensitive to fluctuations. With (O) := Tr(Op), cr 2 yz = 1, and (a x ) = 2Re(pi2), 
(a y ) = -2Im(pi 2 ), (a z ) = 2p u - 1, we obtain: 

A 2 X := (a 2 x ) - (a x ) 2 = 1 - (2Re(p 12 )) 2 , (64) 

:= (a*) - (* y ) 2 = 1- (2Im(p 12 )) 2 , (65) 

A 2 := (a 2 ) - (a z ) 2 = 4 Pll (l- Pll ) . (66) 

While A 2 ]_ 5 f or ry ^ 0, which would correspond to a completely mixed state in QM, A 2 
can have oscillatory contributions with growing amplitude oc exp Cl±t, see Eqs. (|56p -(|59 p . - Also 
note that the would-be QM energy expectation varies in time, a(a z ) = a(2pn(0) — 1) exp(— "y 2 t), 
approaching the stationary mixed-state value zero. 

To summarize: Modifying a QM two-state model by applying a general linear perturbation, 
however small - which preserves hermiticity and trace normalization of the density matrix - 
opens the possibility of a larger state space, reflected in a doubling of the number of degrees of 
freedom. I 12 l 

We emphasize that the QM evolution becomes exponentially unstable when such perturbation 
is introduced. Generically, this spoils the standard interpretation of the eigenvalues of the density 
matrix as probabilities and can result in interesting oscillatory effects, as we have seen. 

It is an interesting question, whether there is a "classical" formulation of general linear 
dynamics, suggested by a doubled number of certain degrees freedom and "anomalous 
probabilities" that seem invariably to appear - and which remind of similar phenomena 
encountered when relating quantum and classical mechanics, e.g., via the Wigner function. 



5. Conclusions 

In this article we have summarized our earlier derivation of a path integral for classical 
Hamiltonian systems based on an ensemble description [32\ [37] . 

This leads us to point out the correlation properties of classical dynamics in parallel to the 
quantum mechanical ones and to identify characteristic similarities and differences. In particular, 
it seems useful to distinguish intra- and interspace entanglement. 

The former has been held characteristic of quantum mechanics and a feature of superpositions 
of tensor product states. The latter concerns classical mechanics only; it correlates Hilbert space 
states and their duals. However, surprisingly, for anharmonic potentials or interactions, classical 
mechanics additionally shows intraspace entanglement, as in quantum mechanics. 

As a first application, we propose a new formulation of hybrid dynamics, i.e., based on a 
hypothetical direct coupling between quantum and classical objects. This may be of practical 
as well as foundational interest. 

12 I.e., a doubling of the initial conditions for the off-diagonal matrix elements here; it will be interesting to see, 
whether this extends to all degrees of freedom in a more general model. 

13 One might speculate on the relevance for flavour oscillations. 



Finally, we study a generalization of quantum evolution, general linear dynamics, where the 
evolution is generated by a superoperator that preserves hermiticity and trace normalization of 
density matrices. We argue that one cannot "leave QM" without giving up one of the three 
defining properties of density matrices, to be Hermitian, normalized, and positive-semidefinite. 

In the most simple case of a two-state system, we solve such dynamics explicitly. We 
show that the corresponding von Neumann dynamics becomes exponentially unstable under the 
influence of a general linear perturbation. Most interestingly, it leads to the appearance of 
"anomalous probabilities" and an enlargement of the state space, possibly pointing towards a 
sort of prequantum dynamics. - We intend to study more complex objects consisting of such 
two-state systems as building blocks which interact. This may be useful in trying to understand 
quantum mechanics as an emergent phenomenon pQ. 
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